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Abstract. Initial results, techniques, and rationale for a near-infrared survey of evolved emission-line stars toward 
the Galactic Center are presented. We use images taken through narrow-band emission-line and continuum filters 
to select candidates for spectroscopic follow-up. The filters are optimized for the detection of Wolf-Rayet stars 
and other objects which exhibit emission-lines in the 2 /im region. Approximately three square degrees along the 
Galactic plane have been analyzed in seven narrow-filters (four emission-lines and three continuum). Four new 
Wolf-Rayet stars have been found which are the subject of a following paper. 
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1. Introduction 

Optical surveys within our Galaxy are severely hampered 
by dust obscuration; therefore complete samples must be 
obtained with longer wavelength studies. Here we de- 
scribe a survey procedure for evolved massive stars in the 
Galactic plane at K-band wavelengths, where the extinc- 
tion is 10 magnitudes lower than for traditional V-band 
surveys. Our scientific driver is the discovery of young stel- 
lar populations in our Galaxy through the detection of 
evolved massive stars. These stars have strong emission 
lines, which makes them relatively easy to detect using 
narrow-band filters. 

Massive stars drive the evolution of galaxies 
through powerful stellar winds, large mass ejec- 
tions, and explosive deaths. These mechanisms are 
the dominant source of energy input into the ISM; 
thus the structure and composition of the ISM in 
most galaxies is largely determined by the mas- 
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also essential contributors to the chemical evolution of 
their host galaxy, ejecting material enriched in helium, 
carbon, and nitrogen during their lives, and depositing 
elements heavier than nitrogen in their final eruption as 
SNe. 

As a massive star evolves, its spectrum becomes dom- 
inated by emission lines, arising either in a dense stellar 
wind, or in circumstellar material produced by mass loss. 
The presence and strength of individual lines are clues to 
the star's evolutionary state and atmospheric structure. 
Among evolved massive stars with such spectra, emission 
lines are most pronounced in Wolf-Rayet (WR) stars and 
in the Luminous Blue Variable (LBV, or S Dor variable) 
stage, where they shed large (1 to 10 Mq) amounts of 



time (~ 10 4 yrs) ( 
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WR stars have lifetimes < 10 Myr, and thus are excel- 
lent tracers of recent star formation, and so also Galactic 
structure. They are also important in our quest to un- 
derstand how star formation proceeds. For example, most 
of the previously known WRs are relatively isolated or 
in OB associations, but recent searches in the IR have 
found a plethora of these objects in compact clusters near 



the Galactic center (Blum et al. 1995, Krabbe et al. 1995 
Nagata et al. 1995| , [Figer et al. 19994 [Blum et al. 200lf 



2. Life Cycle of a Massive Star 



In the 'Conti scenario' 
is determined by the 



the evolution of a massive star 
amount of mass loss it under- 



goes (Conti 1976, Maeder & Conti 1994). This mass loss 



is driven by the combination of the star's intrinsic lumi- 
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nosity and opacity due to metal lines. The more lumi- 
nous a star is, the greater the outward force, and the 
more metals present in the wind, the greater the opac- 
ity and hence the ability to drive a mass-losing stellar 
wind. Our picture of massive stellar evolution is rapidly 
changing, as new models including rotation are devel- 
oped ( [Maeder fc Meynet 2000[ |Meynet fc Maeder 2000| ). 



Table 1. Filter Description 



It now appears that rotation is second only to mass 
loss rate in its effect on massive stellar evolution, and 
it may even be more important at low metallicities 



(Maeder fc Meynct 2001) 



The current qualitative overview of massive stellar evo- 
lution is as follows. For a star with an initial mass of 
> 30 Mq, mass loss and mixing on the main sequence 
deplete the hydrogen-rich envelope and reveal the equilib- 
rium products of CNO-cycle hydrogen burning, creating 
a WN star. Additional mass loss and mixing reveals the 
products of He-burning, and the star becomes a WC star. 
The star will end its life as a Type lb or Ic supernova. This 
evolution may or may not be punctuated by eruptions and 
episodes of large mass loss where the star is identified as 
an LBV. 



3. The need for NIR observations 

The ideal place to study massive star formation and evo- 
lution is our own Galaxy, where we can resolve objects 
on small linear scales and have some hope of complete 
samples. It would be of great interest to find all the WR 
stars in our galaxy to test stellar evolution theories at high 
metallicity, and to learn about environments of massive 
star formation. Another important benefit is the study 
of Galactic structure as traced by young star-forming re- 
gions. 

Surveys for emission line stars in the Milky Way using 
the narrow band technique have been done before, but in 
the optical, where extinction by intervening dust is highly 
problematic (Shara et al. 199E). However, these surveys 
firmly establish WR stars as key tracers of star formation, 



both in the Milky Way and other galaxies (Conti 1991). 
By moving to the NIR, we can identify these objects in 
regions where heavy reddening renders them undetectable 
in the optical. 

With a simple model for the distribution of Galactic 
WR stars, Shara et al. (1999), demonstrate the difficulty 
in conducting WR searches in the optical. Because of high 
extinction in the plane, the vast majority (~ 90%) of the 
Milky Way's massive stellar population is hidden from 
view by obscuring dust. They also showed how this could 
be overcome by moving to the infrared, indicating that 
the number of WR stars as a function of magnitude peaks 
at K= 13 — 14. It is clear that NIR surveys within the 
plane of our Galaxy are essential to significantly enlarge 
the known population of WR stars. 



Central A (fim) 


FWHM (yum) 


2.032 


0.010 


2.062 


0.010 


2.077 


0.015 


2.142 


0.020 


2.161 


0.022 


2.191 


0.013 


2.248 


0.024 



Table 2. Prominent Lines in Wolf-Rayet K-band Spectra 



Central A ((im) 



Transition 



Type 



2.0587 

2.0705/2.0796/2.0842 
2.1126/2.1137 

2.1038/2.1152/2.1155/2.1156 

2.1632 

2.1652 

2.166 

2.189 

2.2779 

2.3178 

2.3470 



He I 2s — 2p 
C IV 3p - 3d 
He I 4s - 3p 
C III/N III 8 - 
He I 7 - 4 
He II 14 - 8 
Br 7, H I 7 - 4 
He II 10 - 7 
C IV 15 - 12 
C IV 17 - 13 
He II 13 - 8 



WN, WC 
WC 

WN, WC 
WC 

WN, WC 
WN, WC 
WN 

WN, WC 
WC 

wc 

WN 



Data is reproduced from Table 2 of Figer et al. (1997), see 
references therein. Wavelengths listed are vacuum 
wavelengths. 



4. Survey Description 

In Table 1 we present central wavelengths and FWHMs for 



our chosen set of K-band filters (Blum & Damincli 1999a, 



Homcicr et al. 2002). Four filters are centered on the char- 
acteristic stellar wind emission lines of He I 2.06 ^m, C 
IV 2.08 ^m, H I Br 7 2.166 /xm, and He II 2.189 ^m, and 
the additional three continuum filters are at 2.03 /xm, 2.14 
//m, and 2.248 /xm. Thus each line filter measurement has 
a continuum point to the red and blue. It is extremely im- 
portant to have continuum points to the red and the blue 
of each line filter because of the variation in continuum 
slope caused by dust extinction. 

In Table 2 we list the most prominent emission lines in 
WR spectra, and the filter response curves are overplotted 
in Figures |l]a and [l]b for unpublished K— band spectra 
of WN and WC stars kindly provided by P. Eenens. In 
both, the three continuum filters are overplotted as dot- 
dash lines, and the four line filters as dotted lines. This 
illustrates the sensitivity of the 2.17 and 2.19 filters to WN 
stars, whereas the 2.06 and 2.08 filters are most sensitive 
to WC stars. 

4.1. Observations and Reductions 

In 1996 the survey was begun at the 1.5m at Cerro Tololo 
Interamerican Observatory, and continued in 1997, 1998, 
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Fig. 1. A representative sample of WR stars with filter transmission curves overplotted, WN stars left panel and WC 
stars right panel. Dotted curves represent the line filters at 2.06, 2.08, 2.166, and 2.19 /xm. Dot-dashed lines represent 
the continuum filters at 2.03, 2.14, and 2.248 /im. The He II filter transmission curve is difficult to see in the left panel 
due to the strong He II lines in the WN stars. WR spectra in both panels are kindly provided by P. Eenens. 




Galactic Longitude 

Fig. 2. Map of the observations from 1996, 1997, and 1998. The x and y axes are in Galactic coordinates. Lines of 
constant right ascension are overplotted in intervals of 2 minutes as dotted lines, and lines of constant declination are 
overplotted every half degree as solid lines. 

Fig. 3. The same data as in Figure || overplotted on the 90cm radio image of the Galactic Center region presented in 
LaRosa et al. 2000. Major features are labeled, and the small rectangle indicates the approximate position of Figure 
1 
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320 318 316 314 312 

Galactic Longitude 

Fig. 4. Same as Figure |], but for the 2000 data. 

Fig. 5. This is a sample of our data from 1996. The image shown here was taken with the 2.248 /im filter. The center 
coordinates are 17:46:35.7, -28:42:35 J2000. The y dimension spans 5' and the x dimension spans ~ 19'. South is up 
and east is to the left. Large scale variations in stellar density are apparent and caused by intervening dark clouds. 



and 2000. For 1996, 97, and 98, we used the same telescope 
and instrument configuration: the Cerro Tololo Infrared 
Imager (CIRIM) with the f/8 mode, giving us a 5' x 5' 
field of view and 1.16" per pixel. In 2000, we used the Ohio 
State Infrared Imaging Spectrometer in the f/8 mode for 
a 10' x 10' field of view and 1.16" per pixel. Our images 
are taken in 'strips' composed of 35 — 45 images. We offset 
1/3 of a chip in RA (1996-1997 data) or Dec (2000 data), 
keeping the other coordinate constant. Thus each spot on 
the image strip is a composite of three exposures. Images 
have been obtained over ~ 3 degrees of longitude between 
±0.5 near the Galactic Center. This is shown in Figure g, 
and overplotted on a 90cm image of the Galactic Center 
region in Figure ||. In 2000, we moved outward along the 
plane to I = 316, at the edge of Centaurus looking towards 
the Scutum-Crux spiral arm. In this region we obtained 
~ 1 degree of Galactic longitude between ±1.0 Galactic 
Latitude. This is shown in Figure || A sample of our data 
is shown in Figure |H[ with a small rectangle indicating the 
approximate position on the 90cm radio image in Figure 

I 



Data reduction is performed with the CIRRED pack- 
age of routines in IRAfJJ], written specifically for CIRIM 
and OSIRIS^ reductions by R. D. Blum. The reduction 
steps are as follows. First, we trim the flat-field images if 
needed. For each filter there is a set of images taken with 
the lamps on and off. We take the median of each set of 
'on' and 'off' frames to make a single 'on' and 'off image. 
We subtract the 'off' frame from the 'on' frame to make 
a flat image, which is then normalized by the mean. A 
bad pixel mask is made by comparing the dome flats with 
the lights on and the lights off and using the histogram of 
pixel intensities to distinguish good and bad pixels. This 
bad pixel mask is then used to correct the flat field images. 

Next, the survey images are trimmed if needed, linear- 
ity corrected with IRLINCOR, divided by exposure time, 
fixed with the bad pixel mask, and divided by the flat field 

1 IRAF is distributed by the National Optical Astronomy 
Observatories 

2 OSIRIS is a collaborative project between the Ohio State 
University and CTIO. OSIRIS was developed through NSF 
grants AST 90-16112 and AST 92-18449 
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image. Then the entire stack of images within a strip is 
used to make a sky image for each filter. The images are 
median combined using 'minmax' rejection with approxi- 
mately half the images thrown out to reject contributions 
from stars. The resulting sky image is subtracted from 
each of the individual frames, and a constant is added 
back to maintain an appropriate sky level in the reduced 
image. These images are then assembled into a strip ap- 
proximately one degree long by cross-correlation of the 
overlapping regions on each frame. Finally, we derive as- 
trometric solutions for each of our strips by comparing the 
248 filter images to the 2MASS catalog images from IPAC 
and the IRAF task CCMAP. 



4.2. Photometry 

Each frame is analyzed with DoPhot 
flBchcchtcr ct al. 1993 ). DoPhot identifies, classifies, 
and performs photometry on objects in an image. It 
makes successive passes over the image, subtracting the 
objects and searching in the next pass for fainter objects. 
The model parameters are found iteratively in the image 
itself, with reasonable first guesses supplied by the user 
in parameter files. 

4.3. Apparent Magnitude Comparison 

Our program has no strict requirement for calibrated 
photometry since emission-line stars are found through 
continuum independent line indices. Furthermore, much 
of the data was taken in non-photometric conditions. 
However, it is of interest to have an order-of-magnitude 
calibration in order to help assess how deep the images 
typically go. We have done a comparison between a typi- 
cal image section of one of our image strips and the 2MASS 
survey images. We compared our instrumental magnitude 
for the 2.06 /mi filter, which is the least sensitive of all our 
filters due to the combination of transmission efficiency 
and filter width. 2MASS magnitudes were obtained from 
the Infrared Science Archive (IRSA) through the GATOR 
query page. Our comparison shows that an instrumental 
magnitude of —3 corresponds roughly to an apparent K 
magnitude of 12 for this strip. 



4.4. Completeness 

To test our sensitivity to emission line stars, we created 
grids of fake stars with continuum magnitudes ranging 
from -2.0 to -7.0 in instrumental continuum magnitude, 
and emission line magnitudes ranging from 0.1 to 1.0 mag- 
nitudes brighter than the continuum. That is, we created 
grids where 100 stars have continuum (2.03, 2.14, and 
2.248 /mi) magnitudes of -2.0, and line (2.062, 2.077, 2.166, 
2.191 /im) magnitudes of -2.1, 100 with the same contin- 
uum magnitudes, and line magnitudes of -2.2, and so on. 
At each step in continuum magnitude, we added 1000 stars 
with varying line magnitude. 



The psf for each filter image was constructed using 
several bright and relatively well-isolated stars in the im- 
age and the IRAF tasks 'pstselect' and 'psf. Images were 
then created using 'mkobjects'. These images were ana- 
lyzed in exactly the same way as a real image. In Figure |^ 
we plot the fraction recovered at increasing brightness of 
the continuum magnitude. There is a dependence on over- 
all brightness, but no significant dependence on emission 
line strength. This is similar to what one would expect for 
a broad-band color. 

The survey is limited by crowding (especially in the 
GC region) and the relatively coarse angular resolution 
element provided by the 1.5m telescope. Given the order- 
of-magnitude comparison made above with the compar- 
ison to 2MASS, the 2.06 /im filter the data are > 75 % 
complete at K = 12 mag. However, the photometric accu- 
racy decreases at faint magnitudes. This is shown in Figure 
^ for the same stars that are plotted in Figure [| The line 
shown here is the 2.062/im line, but all lines (2.077, 2.166, 
and 2.191 /im) show very similar behavior. 

4.5. Candidate Selection 

In this section we detail our candidate selection procedure. 
As a first step, the object coordinates must be transformed 
to a single system. We have arbitrarily chosen to use the 
longest wavelength filter image as our template. Thus, we 
scale and offset the coordinates for objects in the other 
lists to the 2.248 /im object list. This is accomplished 
with the 'transform' package which accompanies DoPhot. 
Transform uses a triangular search routine to match stars 
in different images by comparing the relative scales and 
orientations of stars in groups of three. Hereafter, we refer 
to the task as 'offset'. 

We look at each emission-line filter separately, mean- 
ing we consider sets of three filters each time, the line 
and its two continuum filters (2.03 and 2.14 for 2.06 and 
2.08; 2.14 and 2.248 for 2.17 and 2.19). The raw output 
of DoPhot contains 10,000 to 30,000 objects per image, 
and this output must first be gleaned for stars with good 
statistics. The first cut is on object type. In the DoPhot 
output each object is assigned a number indicating the 
whether it is a star or a galaxy, and how well it thinks 
it can determine the magnitude. There are three object 
types we are interested in: l's are single stars with good 
statistics, 7's are single stars with decent statistics, and 
3's are fit as members of a blend of 2 single stars. During 
the transformation to 2.248 /im coordinates, 'offset' keeps 
track of the object types in both lists, and records it in 
the output file. We keep only combinations of single stars 
(such as ll's, 17's, 71's, or 77's) and doubles (33's). We 
throw out stars which have been classified as a member of 
a double in one image and a single star in another, as the 
photometry in this case is very unreliable. 

After cleaning the lists for object type, we perform a 
cut on photometric error. The error is calculated simply as 
the errors of the individual magnitudes added in quadra- 
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Fig. 6. The fraction of stars recovered in artificial star experiments for the four line images and its dependence on the 
continuum magnitude of the star. See text. 



ture. The error cut is a free parameter in our selection 
routine, but we usually select it to be between 0.1 and 0.3 
magnitudes. 

We now have three lists of objects, cleaned on the basis 
of object type and error. The coordinates in these three 
lists are now matched, and a 'line index' is computed as 
the magnitude difference between the measured value and 
the expected value derived from a linear interpolation us- 
ing the two continuum filters magnitudes. 

Essentially all stars in the images should have a con- 
stant line index (which is not zero due to filter transmis- 
sion differences and sky transmission variations), that we 
can subtract to set the mean line index to zero. When the 
nights were not photometric, systematic line index varia- 
tions as a function of position (time) along the strip occur. 
These are relatively straightforward to account for since 
they produce large scale changes in the line indices with 
position on the image. In our analysis we subtract these 
variations by calculating a "mean neighbor index" for each 
star, using a bin centered on the star, with a size between 
50 — 100 pixels. Smaller bins produce output with less scat- 
ter, but the bin must be large enough to contain stars to 



compute a significant "neighbor index" . This mean value 
is then subtracted from the individual objects. 

We also consider the scatter in line index as a function 
of the short wavelength continuum magnitude. An exam- 
ple is shown in Figure |8|, where the standard deviation in 
line index is fit as a linear function of short wavelength 
continuum magnitude. Our final selection is done by con- 
sidering the line index, the error on this index, and the 
scatter in line index at this magnitude. We employ a cri- 
terion similar to the 'S' parameter employed by Damineli 
et al. (1997). 

An example of our photometry and candidate selection 
is shown in Figures || and [HJ In these figures, candidates 
were chosen as those with k > 0.5, and these candidate 
emission line objects are overplotted as asterisks. One is a 
confirmed WR star of late WC subtype (line index ~ —0.6, 
inst. mag. ~ —3.7) (Homeier et al. 2002). In right hand 
panel of each figure is the the same data plotted, but with- 
out correcting for systematic positional variations in line 
index as described above. One can see that the apparent 
scatter is larger for the objects in the right-hand panels, 
and in the case of the 2.08 filter, it affects which objects 
are chosen as emission-line candidates. In this particular 
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Fig. 7. Standard deviation of the difference between the input color index and the measured line index for the 2.06 
line artificial star experiment. The continuum magnitude is shown on the top of each small panel. The error in line 
index docs not depend on line strength, but does depend on continuum magnitude. 

case, a bona-fide WR star would be missed at this cut on 4.6. Image Subtraction 
k. 
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Fig. 8. The scatter as a function of 2.03/zm magnitude is 
fit with a linear function to aid in candidate selection. 



We also use image subtraction as a way to select high 
priority candidates. This step is only performed for those 
objects that are selected as photometric candidates. All 
procedures discussed here are accomplished with standard 
IRAF routines. An image section is selected within +/— 
125 pixels of a photometric candidate's x position, and a 
y pixel range from 1 to 250 pixels. First the images arc 
cross-correlated and shifted to a common position. The 
two continuum filters are averaged and this result is nor- 
malized by the mean. The line filter is also normalized by 
the mean, and then divided by the averaged, normalized 
continuum image. The result is displayed on the screen 
and examined visually. We have found that this provides 
a useful and necessary check on the photometric selec- 
tion. Effects such as bad pixels and residual images are 
easily seen with this method. We will discuss candidate 
selection in further detail in the upcoming paper on the 
spectroscopic follow-up (Homeier et al. 2003, in prep). 
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Fig. 9. 2.06 Filter: An example of our photometric data and candidate selection. The sta rs marked with aster isks are 
emission line candidates, and one is a confirmed WR star (cindex ~ —0.6, mag ~ —3.7; Homeier et al. 2002 ). In the 
right panel is the same but without correcting for positional line index variations; see text. 
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Fig. 10. Same as Figure O, but for the 2.08^m filter. The confirmed WR star is located at index ~ —0.6, mag ~ —3.7. 



5. Estimating the Number of Candidates 
Expected 

5.1. A Smooth Model 

Following the simple model of Shara et al. (1999), we can 
estimate the number of candidates we expect to find in 
the various regions of our survey. We assume a thin, ex- 
ponential disk, where the stellar density follows a radial 
exponential law: 



TV* = N e- {R - Ro)/aR , 

and an exponential dust distribution: 



-{R-R )/a R 



(1) 



(2) 



where R is the Galactocentric distance, Ro is the so- 
lar Galactocentric distance of 8.0 kpc, an = 3.0 kpc from 
Kent et al. (1991). We take & K ,o to be 0.00008 mag pc" 1 
for consistency with extinction measurements near the 



Galactic Center (Catchpolc et al. 199C). This model pre- 
dicts a maximum extinction at the Galactic Center of 
A K ~ 3.4. 



The WR progenitors, the O stars, should follow this 
exponential law, but we must also consider the metallicity 



dependence for WR formation (Maeder 
which modifies cxr. 

WR/O = O.UZ/Za 



Meynet 1994 



(3) 



For young stars and emission nebulae, the metallicity 
variation within the Galactic disk follows an exponential 
drop-off. Smartt & Rolleston (1997) found an oxygen gra- 
dient in the Galactic disk of: 



12 + log(0/H) = -0.07 x R g (kpc) + 9A(dex) 



(4) 



This leads us to the final equation for WR stellar den- 
sity as a function of longitude and distance from the galac- 
tic center: 



TV* = N 0i w Re 



-{R—Ra)/awR 



(5) 



where olwr has the value 2022 pc and No,WR = 
2.2 kpc~ 2 from studies in the local solar neighbour- 
hood QArmandroff fc Massey 1991| , ponti fc Vacca 1990D , 
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Fig. 11. The number of WR stars expected from our sim- 
ple model per degree of longitude as a function of Galactic 
Longitude. 



Massey 199f| ). This relation predicts a total of ~ 2900 WR 
stars in the Galaxy (compared with ~ 2500 in Shara et 
al. 1999). 

We now have an expression for WR density and dust 
extinction, the next ingredient is absolute K magnitudes 
for WR stars, which are not well-known. From previ- 
ous spectral analyses, WC absolute K magnitudes range 
from -3.8 to -5.7 (ftmartt et al. 200l], |Dessart et al. 200C, 



Hillier fe Miller 1999| |Crowther et al. 2002| ). For late- type 
WCs, which are predominant at high metallicity, ab- 
solute K magnitudes range from -3.8 to -5.0. The ab- 
solute K magnitudes of WR134 and WR136, both 
of WN4— 5s type, are -5.47 and -5.64, respectively 
flCrowthcr fc Smith 1996| ). The Galactic star WR105 



(WN9h) has an absolute magnitude of -6.0 derived from 
its apparent K magnitude , its distance and extinction at 
K (|Churchwell et al. 1992|) . 

Given that WN and WC stars come from a range of ini- 
tial masses, they should also have a range in absolute mag- 
nitude. Additionally, the luminosity of the WNL phase is 
greater than that of the WC phase while at comparable 
temperatures. As expected, it appears from observations 
that WN and WC stars have a range in absolute mag- 
nitude, and the WN star population extends to brighter 
magnitudes. Let us assume a number distribution accord- 
ing to a Salpeter law, that WN stars are distributed be- 
tween K —4.0 and —7.0, WC stars between —3.0 and —6.0, 
and that the WC/WN ratio is 1.2 in the inner galaxy (from 
the em pirical relation for WC/W N and O/H in the Local 
Group, Massey fc Johnson 1998 ). 



5.2. Smooth Model Prediction 

With this model, we can make an estimate of the number 
of candidates we can expect in our data set. In Figure ^ 
we show the number of candidates expected per degree of 
Galactic longitude as a function of longitude, for varying 
apparent magnitude limits. 



Our images cover ~ 3 degrees of longitude, but ~ 15% 
of that is lost by being on the edges of our strips, and 3 
strips are of very low signal-to-noise due to clouds. This 
gives us a total of 2.3 degrees of longitude, covered between 
±0.5 galactic latitude. 

Within that region our model predicts 18 WR stars 
brighter than K = 11, and 34 brighter than K = 12, or 16 
additional stars. Our completeness tests indicate that we 
detect 95% of emission line stars brighter than K= 11, and 
75% brighter than K= 12. According to this, we should 
detect a total of 28 WR stars in this region. In our second 
set of data, taken in 2000 at I = 316, we cover about one 
degree of longitude between b ± 1.0, and expect 4 WR 
stars. 

So far, we have observed all our high priority candi- 
dates in the GC region. We have not yet observed any 
candidates in the Centaurus-Scutum region. In the GC 
region, we have found four new WR stars. These four 
stars should be added to those already found in the region 
through pointed observations. The majority are found in 
two young clusters, the Arches cluster (Nagata et al. 1995; 
Cotera et al. 1996), and the Quintuplet cluster (Figer et al. 
1999), though some stars have been found outside clusters 
as in the present case (Cotera et al. 1999). A full discus- 
sion of our detection efheieny and description of the newly 
discovered WR stars will be given in our upcoming paper 
on spectroscopic follow-up observations. 

Preliminary spectroscopic results are given by 



Homeier et al. 2002 



6. Other Emission-line Objects 

Along with WR stars, we may detect other emission line 
objects in our survey. As mentioned previously, LBVs 
are related to WRs, and also have strong emission lines 
in the infrared. Be and B[e] stars have Br7 lines, but 
much weaker than for WRs or LBVs. Therefore only the 
strongest emitters may be detected. 

Another possibility is WR central stars of planetary 
nebulae. The stars themselves have very weak lines in this 
region, but their surrounding nebulae have lines of He I 
at 2.06 fj,m and B7 at 2.166 /im. 

We can also detect very young O and B stars still en- 
shrouded in their natal gas and dust. These objects form 
compact, or ultra-compact H II regions, and later, when 
they are more revealed, can exhibit emission from cir- 
cumstellar disks. They will be visible in nebular emission 
lines of Br 7 and He I2.06/im (Blum fc Damincli 1999b, 



Hanson et al. 2002). Along with the WR stars which we 



have already detected (see above), we have identified sev- 
eral compact Br7 sources. These will be discussed in our 
following paper detailing the spectroscopic observations of 
our candidate objects. 

7. Summary 

We have completed a near-infrared survey of approxi- 
mately three square degrees toward the inner Galaxy in- 
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eluding a larger region centered on the Galactic center 
(GC) and a smaller region towards the Centaurus-Scutum 
arm. Our survey uses four line emission filters and three 
continuum filters in the K— band. The need for a con- 
tinuum filter on both the red and blue side of each line 
is driven by the large and variable extinction toward the 
inner Galaxy. 

We have recently completed spectroscopic follow-up 
for imaging data taken in 1996-1998 centered on the GC. 
We discovered four new WR stars. The complete results, 
including spectra, coordinates, and finding charts will 
be published in an upcoming paper (Homeier et al. 2003, 
in prep). Preliminary spectroscopic results are given by 
Homeier et al. 2002. 
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